A LIE THEORETIC INTERPRETATION OF MULTIVARIATE 
HYPERGEOMETRIC POLYNOMIALS 



PLAMEN ILIEV* 



Abstract. In 1971 Griffiths used a generating function to define polynomials 
in d variables orthogonal with respect to the multinomial distribution. The 
polynomials possess a duality between the discrete variables and the degree 
indices. In 2004 Mizukawa and Tanaka related these polynomials to character 
algebras and the Gelfand hypergeometric series. Using this approach they 
clarified the duality and obtained a new proof of the orthogonality. In the 
present paper, we interpret these polynomials within the context of the Lie 
algebra sl^+iCC). Our approach yields yet another proof of the orthogonality. 
It also shows that the polynomials satisfy d independent recurrence relations 
each involving + d + 1 terms. This combined with the duality establishes 
their bispectrality. We illustrate our results with several explicit examples. 

1. Introduction 

We start by introducing some basic notations which wih be used throughout 
the paper. The first important object is the set ICd which parametrizes families of 
multivariate Krawtchouk polynomials. In order to motivate the definition of ]Cd, let 
us first consider positive real numbers {pj}'j^Q such that po+pi + - ■ ■+Pd = 1 and let 
P = diag(po,Pi, ■ • ■ ,Pd) be the corresponding diagonal matrix. Let wo,wi, . . . ,Wd 
denote mutually orthogonal vectors in R"'"'"^ with respect to the inner product 
(a, /3) = a^-PP such that wq = (1, 1, . . . , 1)* and the 0-th coordinate of Wj is 1 for 
j = 1, 2, . . . , d. If we denote by U the matrix with columns wq,wi^ . . . ^Wd then 
Q — U^PU is a nonsingular diagonal matrix whose (0, 0) entry is equal to 1. If we 
set P — pqQ~^, then P — diag(po,Pi,P2, ■ • ■ ,Pd) is also diagonal and 

—PUPU' = Id+i, 
Po 

where Id+i denotes the identity (d+l) x 1) matrix. The above construction is 
closely related to the one by Griffiths [6 . It shows how to build matrices U and P 
starting from P. In the present paper we will continue to work with the matrices 
P, P, U but we would like to put P and P on the same footing. Moreover, we would 
like to work over C instead of R. This leads to the following formal definition which 
extracts the important properties of the matrices P, P and U needed in the paper. 

Definition 1.1. Let ICd denote the set of 4-tuples {ly, P, P, U), where ly is a nonzero 
complex number and P, P, U are (d+l) x (d + l) matrices with complex entries 
satisfying the following conditions: 
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(i) P = diag(po,Pi, ■ • • and P = diag(po,i5i, ■ . ■ ,Pd) are diagonal and po 



(ii) U — {ui_j)o<cij<d is such that uo.j = Uj,o — 1 for ah j = 0, 1, . . . , d, i.e. 

/I 1 1 ... 1 \ 

?7= . ; (1.1) 



\1 UdS Ud^2 ■ • • Ud,dJ 

(iii) The following matrix equation holds 

i^PUPU* - /rf+i. (1.2) 

We note that the points in ICd arise naturally from the so called character algebras 
[T]. In particular, Hecke algebras of Gelfand pairs, or more generally, the Bose- 
Mesner algebras of commutative association schemes are examples of character 
algebras. From the above definition it is easy to see that pj and pj are nonzero 
numbers, such that 

d d 

i=o j=o 

Griffiths |B] used a generating function to construct d- variable Krawtchouk poly- 
nomials for every point k G /Cd and positive integer iV. His construction is as 
follows. For m — (mi, m2, . . . , m^) G Nq and rh = {mi,rh2, . . . ,nid) G Np such 
that mi -|- m2 + • • • + m^ < N, mi + m2 + • • • + ifid < N define polynomi- 
als V{m,rh; K, N) = T'{m,rh) in the variables fhi, . . . ^rfid with degree indices 
mi, . . . ,md by 

1 1 rVim, fh)z'^^ ■ ■ ■ z^""* , 

mo!mi! • • -md! ^ ' ^ ^ ' 

m-i-\ \-md<N 

where mo — N — mi — m2 — • • • — m^ and mo = N — ifii — m2 — • ■ • — rhd. 

Mizukawa and Tanaka [13) gave an explicit formula for V{m, rh) in terms of the 
Gelfand hypergeometric series 

^^"'"^ = ^ i^N)^. , n 

A={ai,j)eMi,N 'l^i,j=i°-i.i i,j=l '-J 

(1.3) 

where uji_j — 1 — Ui,j. In the above formula Aid.N denotes the set of all c? x c? 
matrices A — (dij) with nonnegative integer entries, such that X]f j=i ^* j — 
One advantage of equation ()1.3|) is the transparent symmetry between m and rh. 

In the present paper, we define two Cartan subalgebras for s[d+i(C), denoted 
H and H, which depend on k. We display an antiautomorphism a of s[c;+i(C) 
that fixes each element of H and each element of H. We consider a certain finite- 
dimensional irreducible s[(i+i(C)-module V consisting of homogeneous polynomials 
in (d+l) variables of total degree N. We define a nondegenerate symmetric bilinear 
form {,) onV such that {P.^,r]) = o(/3).77) for ah /3 G sU+iiC) and ^,77 G V. 
We define also two bases for V; one diagonalizes H and the other diagonalizes 
H. Both bases are orthogonal with respect to (, ). We show that when ( , ) is 
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applied to a vector in each basis, the resuh is a trivial factor times one of the 
polynomials 'P(to,to) defined in ()1.3p . Thus the transition matrices between the 
bases are described by the polynomials (jl.3p . From these results we recover the 
previously known orthogonality relation proved in [6l 113) . 

Our approach naturally leads to two sets of d recurrence relations for the poly- 
nomials (|1.3p . parametrized by the Cartan algebras H and H. Equivalently, we 
can interpret these recurrence relations as two commutative algebras of difference 
operators, each generated by d algebraically independent operators diagonalized 
by the polynomials V{m,rh). One of the algebras consists of difference operators 
acting on the variables m and the other one on the variables m. Following the 
terminology established in the literature starting with [3] we can say that the poly- 
nomials V{m, in) solve a discrete-discrete bispectral problem. The two algebras can 
be connected via a natural (bispectral) involution b on Kd defined as follows: 

h:K = {v,P, P, U) b{n) = {v, P, P, [/*). 

This map exchanges the roles m and rh and corresponds to a symmetry between 
the degree parameters and the variables of the polynomials. 

At the end of the paper, we illustrate how our theory applies to multivariate 
Krawtchouk polynomials, which have appeared in different applications in the lit- 
erature. The first example explains how the present work extends our joint paper 
with Terwilliger [TU] related to the bivariate polynomials defined in P|. The next 
example concerns the polynomials defined by Milch more than 40 years ago 
within the context of multi-dimensional growth birth and death processes. We fix 
the matrix P and we exhibit explicit matrices P and U, satisfying the conditions 
in Definition 11.11 Our constructions lead to the bispectral commutative algebras of 
difference operators obtained recently in [5l Section 5.4]. The bispectral involution 
b above corresponds to the mapping f in [5l page 450, eq. (5.22)]. In the last 
example we explain how the polynomials discussed in [2] fit within our framework. 

The results of the present paper yield a family of solutions to Problem 7.1 in 
[TU] . Each of these solutions can be viewed as a rank d generalization of Leonard 
pair. The Leonard pairs are defined and classified in |14| . For more information on 
Leonard pairs see [15]. 



2. Cartan subalgebras of s[d+i(C) 

For i,j e {0, 1, . . . let Cij denote the (d+l) x {d+ 1) matrix that has 
entry 1 and all other entries 0. We denote by H the standard Cartan subalgebra 
of sld+i{C) consisting of all diagonal matrices with basis {4>i, (f>2, ■ • ■ , 0d} where 

(pi = ei,i - -j-^Id+i, for i 1, 2, . . . , d. (2.1) 

Let R denote the (d-|-l) x (d+l) matrix given by 

R = ePu\ (2.2) 

where S C is such that det(i?) = 1. We consider the automorphism Adji on 
s[d+i(C) defined by AdB.W) = R^R'^ for every /3 e sld+i(C). We denote by 0fc 
and Cij the images of 0^ and e^j under Ad^, i.e. we set 

0fc = R(j)kR'~\ = Re^,jR-^ (2.3) 
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for fc = 1, 2, . . . , d and i,j G {0, 1, ... , d}. Let 

7? = span{0i,^2, • • • jf^d} (2.4) 
denote the conjugated Cartan subalgebra of sld+i{C). From (|1.2p and (|2.2p we see 



that 



= 6'PC/, where 



(2.5) 



Using the above equations we can easily expand 0^ in terms of the basis {4>j, e^j }, 
j & {l,2,...,d}, k^l e {0,l,...,d} of s[d+i(C) as fohows 

d 

4>i = v ^ PiPkUiMUtjCkj +^Pi{vpjufj ~ (2.6) 

0<k^l<d 3 = 1 

for i — l,2,...,(i. Likewise, we can also expand (pi in terms of the dual basis 
{h^^kA.j e {l,2,...,d}, fc^Z e {0,l,...,d} ofsld+i(C): 

d 

(f>i^v ^ PiPkUk,iUi^iek,i + '^Pi{vpjV?ji - l)4>j, (2.7) 

0<k^l<d j=l 



for i = 1, 2, . . . , d. 



Next we define an antiautomorphism a on s[(;+i(C) by 

a(^) = P/3*F-i for every /3 £ sld+i(C). 



(2., 



Note that a is an involution, i.e. a o a = Id. Using that fact that P and P are 
diagonal matrices and (II. 2p one can check the following lemma. 



Lemma 2.1. We have 
a(0i) = (l>i 



a(0i) = (pi 



Pi Pi 



for alii = 1,2,..., d (2.9) 
/or allO<i^j<d. (2.10) 



/n particular, a preserves the Cartan subalgebras H and H . 

Lemma 2.2. The Cartan subalgebras H and H together generate s[rf+i(C). 

Proof. Consider 

d ^ 



and 



Using equations (12.21) , (|2.5I) and the explicit form of the matrices P, P and U we 
see that 

00 = i^PU*eo,oPU - -T^Id+i 
d + I 

fPo Pq Po ■■■ Po\ 
pi pi pi ... pi 



d+1 



Id+l- 



\Pd Pd Pd ■■■ PdJ 
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From the last equation it follows easily that for < i 7^ j < d we have 

_ [0j>o]]] " [(f>i, 0o]] 

completing the proof. □ 

3. An s[rf+i(C)-MODULE 

Let xq,xi, . . . ,Xii denote mutually commuting variables. We set x = {xq, xi, . . . , x^) 
and we denote by vo,vi, . . . ,Vd the standard basis of C*^^. Throughout the paper 
we use multi-index notation. For instance, C[x] stands for the C-algebra consisting 
of all polynomials in xo,xi, . . . ,Xd and if A = (Aq, Ai, . . . , Xd) £ Ng^^ then 

x^=x^«a;^•••a;^^ A! = AqIAi! • • • Ad!, 

A| = Aq + Ai + ■ • ■ + Xd, etc. We shall denote by p = (po,Pi, ■ • ■ ,Pd) and p = 
{po,Pi, . . . ,pd) the vectors whose components are the entries of the diagonal matri- 
ces P and P in Definition ll.il 

Consider the representation p : slc;+i(C) qI{C[x]) oi sld+i{C) defined by 

for i ^ j. Thus every element of sld+i(C) acts as a derivation. Let 

1= {AeN;5+i ; |A| ^N} 

and let V denote the subspace of C[x] consisting of homogeneous polynomials of 
total degree N, i.e. 

V = span{x^ : A e I}. (3.2) 

Since the operators in p.ip preserve the total degree of the polynomials, we can 
consider V as an s[d+i(C)-submodule of C[x] and it is easy to see that this module 
is irreducible. In particular, from Lemma 12.21 it follows that there is no proper 
subspace W oi V such that HW C W and HW C W. 
Since 

(j),.x^ = (^X, - -j^^ x^ for i = 1,2, . . . ,d, (3.3) 

dj.x^ = Aja;^+^--^^ for < i 7^ j < d, (3.4) 

we see that 

V = ^Vx, where Vx = spanjx^} (3.5) 

Aei 

is the weight decomposition of V with respect to H. To describe the weight de- 
composition of V with respect to H we need to make a change of variables. Let us 
define x = (io, ii, ■ . . , id) by 

X = xR, (3-6) 

where R is the matrix given in (|2.2p . From equations (|2.3p . p.3p and (13. 4p it follows 
that 

4>,.x^ ^ (^X, ^ x^ for^ = l,2,...,d, (3.7) 

eij.i^ = Aji^+^--^^ for < i 7^ J < d. (3.8) 
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Thus 

y = Va, where = span{i^} (3.9) 
is the weight decomposition of V with respect to H. 

Remark 3.1. We comment on how H and H act on the weight spaces of the other 
one. A pair of elements A and /i in I will be called adjacent whenever A — ^ is a 
permutation of (1, —1, 0, 0, . . . , 0) G C'^^. Then H and H act on the weight spaces 
of the other one as follows. For all A G I, 

mGI Men 

1.1 adj A (1 adj A 

We conclude the section by writing explicit formulas for the entries of x, which 
will be needed later. Using p.6|) . the notations in Definition 11.11 and p.2p we obtain 

d 

io = (3.10) 

0=0 

d 

Xk = Opoxo + d^^UkjPjXj for k ^ 1,2, . . . ,d. (3.11) 

Recall that uJij = 1 — Uij and therefore, subtracting p.lOp from p. Ill) we can 
rewrite equations (|3.11l) as follows 

d 

Xk — xq — O'^^PjUJk^jXj, ior k — 1,2, . . . ,d. (3-12) 



4. A BILINEAR FORM 

We define a symmetric bilinear form (, ) on y by 

= (5„„j— 6*^ for all n,m €l. (4.1) 

Using the explicit formulas for the action of a on the basis Ck^i), j G {1, 2, . . . , d}, 
k I G {0,1, . ■ . , d} oi s[d+i(C) in Lemma I^Tl it is easy to check that 

(/3.e,?7) = (e,a(/3).r7) forall/3Gs[rf+i(C), ^,veV. (4.2) 

We prove next that the vectors i" are also mutually orthogonal with respect to (, ). 

Lemma 4.1. For n,™ G I we have 

=5„,™ — 0^. (4.3) 

Proof. For n ^ m, equation (|4.3p follows immediately from (|4.2p . the fact that a 
preserves H and (|3.7p . Thus it remains to consider the case n = m, i.e. we need to 
show that 
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We prove that ()4.4p is true by induction on ni+n2 + ■■ ■+nd- li ni+n2 + - ■ ■+nd = 
then n = (TV, 0,0, ... ,0) and we need to show that = NW^ /p^ . Taking the 

A^th power of equation (j3.10|) and using (|4.ip we find: 



q2N 



q2N 





2 








> — rP a; 







qG 



(a!)^ 



qGI 



(iV 



1^2 



-p"d 



N 



Po 



completing the proof in this case. Suppose now ni+n2 + - ■ ■ + nd > 0. Then rij > 
for some j G { 1 , . . . , d} and we consider the equahty 

From ((XSI it fohows that the left-hand side equals | | p. From (pHI)) 

and p.Sp we see that the right-hand side is ^("-o + Therefore, we have 



<~nu2_P0 nj 



Pj ("0 + 1) 



\~n+vo~Vj 1 12 



completing the proof by induction. 



□ 



5. The Krawtchouk polynomials and s[d+i(C) 

We are now ready to prove that the transition matrices between the bases {x" : 
n G 1} and {i" : n G 1} can be expressed in terms of the multivariate Krawtchouk 
polynomials ()1.3p . To state the theorem we shall use the following convention: for 
a vector w = {wo,wi, . . . , Wd) G C''+"'^ we denote by w' — {wi, . . . , Wd) G C^. 

Theorem 5.1. With the above convention we have 

^ = my^V(n',ii')^x'' for fie I (5.1a) 



— ^ N\Y,V{n' ,fi')^x^ for nel. (5.1b) 



Proof. Let us denote uji — {u!i^i,u!i,2, . ■ . ,^i,d). Using equation p. 121) we find 
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If we set a — (oi,!, ai,2, • ■ • , di^d) and denote by A the matrix with entries (aij), 
then multiplying the above equations we obtain 

0N 



AeMd,N ^ ' j=i j=i 



d a-i, 

n 



(5.2) 



X 



Using now p.lOp we get 



E ^ — ;r' ' n(MP. (5.3) 

l7l=W-E-,, = ia.,i 

Note that 

If we change the variables as follows no — 70 and rij = 7^ + X^iLi J — 

1, 2, . . . , d and we substitute the last formula in (|5.3p we obtain 



/ ^-Ef,j=i Qi.j 2 iV! 

The proof of (|5.1ap now follows immediately by plugging (|5.4p in (|5.2p . The proof 
of (jS.lbp can be obtained along the same lines, or one can use the duality. □ 

Corollary 5.2. For n,n G I we have 

P{n',h')^-^{x\x"). (5.5) 



As an immediate consequence of Theorem 15.11 and Corollary 15. 2[ we obtain a 
new proof of the orthogonality relations established in [6l [13] . 

Corollary 5.3. For h, k El we have 

N\Y^V{n\n')r{n\k'f-=5,;,j^^^ (5.6a) 



and similarly for k £ I we have 



Nlj:nn\n')nk',n'f-=S^^,j^^. (5.6b) 

Proof. Equation (|5.6ap follows if we consider the inner product of each side of 
equation ()5.1ap with x'' and then we use (|4.3p and (j5.5p . The proof of (|5.6bp is 
similar. □ 
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6. BiSPECTRAL DIFFERENCE EQUATIONS 

III this section we show that the Cartan subalgebras H and H lead naturaUy to 
recurrence relations for the multivariate Krawtchouk polynomials. We work below 
with and we denote by wi, W2, • • • , "I'd its standard basis. 

Theorem 6.1. Fix m = (mi, TO2, . . . , mj) & Nq and rh — (mi, 7712, . . . , rhd) e Ng 
such that |m| < N and \rh\ < N. Then for every i G {1, 2, . . . , d} we have 



N \ 



m — vi) 



d 

d 

uy^^PiPkUk.ij^ - |m|)7'(m,m + Vk) 

k=l 

^pr{vpju\^ - 1) im.j - j V{m,m) 

j=i 

V ^ piPkUk,iUi^^miV{m,fh + Vk~vi), (6.1a) 

l<k^l<d 

N \ 

^^i j P(m,m) = 'Y^piUi^imiVim - vi,m) 

d 

vy^^PiPkUi,k{N - \m\)P{m + Vk, rh) 



+ '^P^i^P]ulJ - 1) {ruj - j^j V{m,m) 
j=i \ 

+ I' ^ piPkU.LkUi,im'P{m + Vk - vi,'ih). (6.1b) 

l<k^l<d 

Proof. For m, ifi E Nq such that |m| < N and |m| < N, we consider n = {N — 
|m|, mi, m2, . . . , m^) G I and n = {N — |m|, mi, 7712, . . . , fhd) G I. From Lemma |2. II 
and it follows that 

(0,.x",i") = (x",0,.i"). (6.2) 

From (13. 3p and Corollarv l5.2l we see that the left-hand side of (|6.2p is equal to N\ 
times the left-hand side of (|6.1ap . Now we use (12. 7p and Lemma 12.11 to evaluate 
the right-hand side of equation (|6.2p . which gives N\ times the right-hand side 
of ((5?Ta|) . The proof of (|6.1b|) is similar. □ 

Remark 6.2. Theorem l6 . 1 [ establishes the bispectrality of the polynomials V{m, rh). 
Indeed, equations (j6.1ap show that the polynomials V{m,rh) are common eigen- 
functions of d difference operators acting on the variables mi, m2, . . . , m^, with 
coefficients independent of mi, m2, . . . , m^, while equations (I6.1bp indicate that 
the same polynomials are common eigenfunctions of d difference operators acting 
on the variables mi, m2, . . . , m^, with coefficients independent of mi, m2, ■ . ■ , rhd. 

Formulas (|6.ip show that the polynomials Vim^fh) form a remarkable basis. 
A random basis of polynomials orthogonal with respect to the the multinomial 
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distribution will not satisfy such simple recurrence relations, see [4] for the general 
theory. If we add equations ()6.1ap for i — 1,2, ... ,d and use the matrix equation 
(|1.2[) we obtain 

d d 
— \m\'P{m, m) —po fhiV{m, rfi — vi) + ^^pi(iV — |m|)7^(m, m + vi) 
1=1 1=1 

+ ^ PkfhiV{m,m + Vk-vi) 

i<k^i<d ^ ' ' 

( ^ \ 

+ I "f^jPj PoiN — |m|) — N V(m,fh). 
Vj=i j 

Note that the coefficients on the right-hand side in the above partial difference 
equation depend only on p, N and m, while the eigenvalue on the left-hand side 
depends only on the total degree \m\ (if we think of V{m,fh) as a polynomial of 
to). This equation is universal since every basis of polynomials orthogonal with 
respect to the multinomial distribution 

N—\m\ rhi m2 rhd 

Pq Pi P2 ■■■Pd 



wiU satisfy (|6.3|) . see [TT]. 



7.1. Let d = 2 and define 



[N — |to|)!toi!to2! • • • 



7. Explicit examples 



T (Pl +P2)(P1 +P3) , (Pl +P2)(P2 +P4) 

"1,1 = 1 ; ; ; -, ^1,2 = !- 



Pl(pl+p2+p3 + p4)' ■ p2(pl +p2 + P3 +P4) 

, (P1+P3)(P3 + P4) ^ (P2 +P4)(P3 +P4) 

"2,1 = 1 7 T, U2,2 = 1 



P3(P1 +P2 +P3 +P4)' ' P4(pl +P2 + P3 + P4) 



Plp2(pl + P2 + P3 + P4) P3P4(P1 + P2 + P3 + P4) 



(pl + P2)(pl + P3)(P2 + P4) ' (pl + P3)(P2 + P4)(p3 + P4) ' 

plp3(pl -Kp2 +P3 +P4) . P2p4(pl +P2 +P3 +P4) 

Pl = T. , ■ N/. , . N/. — -P2 - 



and 



(pl + P2)(pl + P3)(P3 + P4) ' (pl + P2)(p2 + P4)(p3 + P4) ' 



Po = P0 = 1 - Pl -P2, 



where the numbers {pj"}i<_;"<4 are essentially arbitrary, although certain combina- 
tions are forbidden in order to avoid dividing by 0. A straightforward computation 
shows that (II. 2p holds and we obtain the polynomials defined by Hoare and Rah- 
man The constructions in the present paper correspond to the ones in jlO) . 
The bispectral involution b in this case amounts to exchanging the parameters p2 
and p3. We note also that independently of [TO], Griinbaum and Rahman [5] used 
different methods to derive an interesting recurrence relation for these polynomials 
which was suggested in [7]. 
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7.2. Within the context of orthogonal polynomials one is often interested in ex- 
plicit formulas for the polynomials in terms of the parameters which appear in the 
measure. We construct one such example which gives the multivariate Krawtchouk 
polynomials discovered by Milch [12_. 

Let us fix nonzero complex numbers po,pi, . . . ,pd such that po+pi + ■ • ■+Pd = 1- 
We define pQ,pi, . . . ,pd in terms of po,Pi, ■ • ■ iPrf as follows: 



Pk = ,^ _ ^ forfc = l,2,...,d. (7.1b) 



Po^Po (7.1a) 
PkPo 

(i-E-=iP.)(i-E-;iP,) 

Next we define a (d + 1) x (d + 1) matrix U = {uij) with entries 

Uij = (5o,j, when i < j, (7.2a) 

Uij = 1, when i > j, ('''.2b) 

uo,o = 1, (7.2c) 

Ui, = ~ 2^k=iPk ^ for t= l,2,...,d. (7.2d) 

Pt 

Thus, J7 is a matrix of the form (jl.ll) where the remaining entries are O's and 
I's above and below the diagonal respectively, and the diagonal entries are given 
in (j7.2dl) . One can check that with the above notations equation (|1.2p holds and 
therefore we obtain explicit formulas for all constructions in the paper. Up to a 
permutation of the variables and the parameters, this choice leads to the multivari- 
ate Krawtchouk polynomials considered in |12j and the corresponding bispectral 
difference operators constructed in [51 Section 5.4]. 

7.3. Fix now q ^ {0, 1} and define pa,Pi, ■ ■ ■ ,Pd as follows: 

Po = q-" (7.3a) 

Pk = q-''+''-\q - 1) for A: = 1, 2, . . . , d. (7.3b) 

We denote 

P = P = diag(po,Pi, ■ • ■ ,Prf), (7.4) 
and we consider also a (d + 1) x (d + 1) matrix U — {ui,j) with entries 

Uij- = 1, when i + j < d, (7.5a) 

when i + j = d + 1, (7.5b) 



Uij = 0, when i + j > d + 1. (7.5c) 

With the above notations, we see that (|1.2p holds leading to the family of multi- 
variate Krawtchouk polynomials introduced in [5] . 
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